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Abstract. - Background charge rearrangements in metallic single-electron transistors are 
modelled in two-level tunnelling systems as a Poisson process with a scale parameter as only 
variable. The model explains the recent observation of asymmetric Coulomb blockade peak 
spacing distributions in metallic single-electron transistors. From the scale parameter we esti- 
mate the average size of the tunnelling systems, their density of states, and the height of their 
energy barrier. We conclude that the observed background charge rearrangements predomi- 
nantly take place in the substrate of the single-electron transistor. 



Introduction. - The metallic single-electron transistor (SET) is a possible building 
block of future electronics, based upon the controlled transfer of individual charges onto and 
off small isolated electrodes, called islands. However, due to the extreme sensitivity of SETs 
to charges close to these islands, a static background charge arrangement is essential for 
proper device operation H[|]. On the other hand, this very high sensitivity makes the SET 
an excellent device to investigate background charge rearrangements, e.g. in nearby two-level 
tunnelling systems (TLTS) |,§. One issue of this paper is to determine the location of TLTSs 
within the device. 

In recent experiments |7|,|^] on Al/A10x/Al SETs it was shown that the distribution of 
nearest-neighbour spacings (NNS) between Coulomb blockade oscillation peaks is influenced 
by such charge rearrangements: they generate a pronounced tail in the NNS distribution 
towards smaller spacings (see left inset of fig. ^). Here, we propose a quantitative model for 
this tail. We argue that background charge rearrangements can be interpreted in terms of a 
Poisson process. This way, we can fit the shape of the NNS distribution tail, with a single 
scale parameter. From this fit, we find an average hopping distance in the TLTS of about 4nm 
and conclude that the charge rearrangements do not occur in the tunnel barriers of the SET. 
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We can furthermore estimate the density of states of TLTSs and the barrier height between 
the two TLTS states. 

Experimental observations and statistical model. - The NNS distribution p(AV g ) in fig. [I] 
is strongly peaked at AV g = e/C g , but displays a remarkable tail towards lower AV g as well, 
which contains up to one third of the total data points. Here, V g and C g denote the gate voltage 
and the gate capacitance, respectively. The spacings observed in the low-AV^, tail are stable: 
repeated and reversed gate voltage sweeps reveal the same statistical distribution of events. 
It has been attributed to the rearrangement of background charges in TLTSs [||, fig. ||a). A 
further important experimental result is that the NNS distribution does not depend on V g : 
it was found Q that the standard deviation of the n-th neighbour Coulomb peak spacings is 
proportional to y/n. 

The clearly observed hysteresis of the Coulomb peak pattern when reversing the gate 
scan direction Q suggests the existence of metastable configurations, which explains the high 
stability and reproducibility of the experimental results. Neither the Coulomb blockade peak 
positions nor the underlying smooth 1// noise spectrum depend on V g or the scan speed. 
Hence, other excitations which could spontaneously flip a TLTS (like thermal activation) are 
much weaker compared to the parameters which we control in the experiment. Note that 
these experiments are different from dynamic and frequency dependent background charge 
noise measurements |], pi PM2] • 

The focus of our model lies on the p(AV g ) distribution generated by TLTS fluctuations. 
Since the AV g distribution is stable, a single sweep of V g already reveals it completely. There- 
fore, our study is restricted to a single gate sweep. During this sweep, the electric field around 
the island electrode retains its direction thus causing all background charges in TLTSs to 
'hop' either towards or away from the island electrode. Therefore, the total amount of charge 
Qi influenced on the island is monotonic as a function of V g and is tested periodically (with 
respect to V g ) by the position of the Coulomb peak. Background charges in TLTSs contribute 
to the variation of AV g by influencing additional charge on the central island. 

We display the situation of the experiment in fig. ^ schematically. While sweeping V g , 
TLTS switching events lead to stochastic jumps of the island charge Qi, see fig. |&). Due 
to the discrete technique of measuring the Coulomb spacings, the Qi jumps are recorded at 
almost regular V g intervals. From a mathematical point of view, we can consider V g as a 
function of Qi just as well. In this picture, see fig. |]c), we observe steps of V g of (almost) 
uniform height AV g which occur at arbitrary values of Qi. 

The sequence of Qi values, at which V g jumps, forms a stochastic process where the 
charge Qi takes over the role of the time in the common notation of stochastic processes. 
Assuming statistic independence of the jumps (see below), the process can be identified as 
Poisson process |l3) . The single requirement for a Poisson process is the independence of the 
individual events. The process is described by a single scale parameter I/7. The probability 
that n jumps occur before the charge Qi is accumulated on the island, is given by [|l3| 

Pn(Qi) = e*p(-7Qi). (1) 

The Poisson process is non-stationary, which in our context corresponds to an average 
number of jumps (n) = yQi, which depends on Qi. jQi also corresponds to the square 
standard deviation, i.e. ((n)) = (n 2 ) — (n) 2 = jQi. 

In the framework of the Poisson process it is possible to describe the tail of the p{AV g ) 
distribution function. For that purpose we transform p( A V g ) intop(AQj), which follows from 



H.-O. MULLER et al: MODELLING BACKGROUND CHARGE REARRANGEMENTS 



3 



AQ l = e-AV g C g . (2) 

p(AQi) is peaked at AQi = and displays a tail into the positive, see right inset in fig. |l|. 
The probability distribution of the charge difference AQ; between two jumps corresponds to 
the waiting time distribution of the conventional Poisson process. Therefore 

p(AQ t ) = - dP f^ l) = 7 exp(- 7 AQ,) (3) 

is the tail of the distribution function p(AQi) with po(AQi) as defined in eq. ([!]). Finally, 
the mean value and square standard deviation of AQi also follow from the distribution function 
p(AQ 4 ): 

1 



(AQi) = V«AQi» = -■ (4) 
7 

Discussion. - It has to be noted that the presented stochastic model is only valid for the 
description of the tail of the p(AQi) distribution, but not its maximum. While the tail is due 
to events arising from randomly distributed TLTSs, for which the presumptions of the Poisson 
process will be proven to hold, the maximum reflects the regularly spaced Coulomb blockade 
peaks and is deterministic in nature. Fortunately, these two ranges can easily be distinguished 
in our experimental data (see fig. |l|). In order to separate both domains, a phenomenological 
cut is made at AQ ijmin = 0.01 e, which corresponds via (||) to AV^ max = 0.99 e/C g . 

Our model provides us with three independent methods of determining the value of its 
sole parameter I/7: (i) a fit to eq. (||), (ii) the average (AQj), and (iii) the deviation ((AQi)} 
according to eq. (Q). For the data presented in fig. [I] we find the respective values: (i) 1/7 = 
(0.046 ± 0.002) e, (ii) 1/7 = 0.057 e, and (iii) l/ 7 = 0.043 e. The discrepancy between the 
average obtained from (ii) and the other values can be shown to be due to the AQi !m i n cutoff. 
However, the correspondence between (i) and (iii) is sufficient to demonstrate that the TLTS 
events are statistically independent because it reflects a unique feature of the Poisson process. 
This process provides an adequate description of p(AQ^) and in turn p(AV g ). Measurements 
on five different SET devices of equal geometry and on the same substrate result in a mean 
1/7= (0.047 ± 0.003) e. 

The value of 1/7 provides a charge scale of the influence of TLTSs on the island charge Qj. 
Thus, it describes the material properties and consequences of the sample manufacture. A 
small value of 1/7 corresponds to a short tail of the distribution functions p(AQi) and p(AV g ), 
which in turn correlates with few TLTS transitions per e/C g interval and, thus, rather stable 
device operation. 1/7 also scales with device dimensions, as will be shown below. 

A simple geometric argument allows us to deduce a mean tunnelling distance 8 of the 
TLTS from 1/7. It is based on the "constant capacitance model" |l4|]: charge rearrangements 
within TLTSs induce potential shifts only, but they do not alter any capacitances. This is 
appropriate here, since only metal electrodes with high electron density and individual charges 
at the TLTS sites are considered. The assumption of constant capacitances is also supported 
by the experimental data, since the peak position of p(AV g ) is independent of V g and the 
bias voltage (see figs. Id) and 3&) in ref. 0, respectively). In addition, screening of the single 
TLTS charges is neglected, which is reasonable considering the low concentration of quasi-free 
electrons around them, i.e. in silicon oxide or aluminum oxide. The treatment is simplified 
by assuming a constant field between gate and island, and by assuming that only one electron 
is rearranged in a TLTS switching event. 
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In this case, Green's reciprocation theorem |15| can be applied to the setup displayed in 
%• la), 

n n 
3=1 3=1 

where {Qj} and {Q'j} are two charge configurations in a system of conductors and {Vj} and 
{Vj} are the respective voltage configurations. 

{Qj} shall be the charge configuration before and {Q'j} that after charge transfer. The 
influence on the gate electrode is negligible, i.e. Qi = Q' x and V\ = V[ = V g . Furthermore, 
we use a single charged TLTS, 

Q 2 =e Q' 2 = Q 3 = Q' 3 = e. 

In our model the charge on the island is conserved, Q4 = Q' 4 = Qi- The island potential, 
however, shows a slight increase, VI = V4 + AV4. This yields with eq. (||) 

m = 4 {V3 ~ vi) - (6) 

For symmetry reasons the mutual contributions to the potentials V3 and V{ cancel and 
(V3 — Vg) is solely determined by V g , 

Vs-Vl = ^V g (7) 

using the gate-island separation D and the TLTS size 5 (see fig. ||a)). In general, eq. (0) 
looks different for a different electrode geometry, but a characteristic length scale correspond- 
ing to D can be found and a similar formula will hold in most cases. The equation gives 
an upper estimate of the voltage difference since it assumes the TLTS to live in the high- 
field region. This region experiences the widest energy range of TLTS transitions and thus 
contributes most events to the statistics. 

Finally, Qi = C g V g together with @ and (Q) yields 

This value is independent of V g , as found in the experiment, due to the cancellation of the 
linear V g dependence of Qi and (V3 — V£). 

Fixed capacitances allow us to relate AV4 of eq. (§1) to AV g of eq. (||) and to establish a 
link between I/7 of (||) and 6 of (||). Given the regular Coulomb peak separation e/C g , the 
link is provided by 

e/Cg e e-f 

Hence, 7 scales linearly with D, i.e. the NNS tail is experimentally observable only for 
sufficiently small D. However, the result (^|) is independent of the distance between the 
TLTS and the island. This is due to the assumption that the charge of the TLTS influences 
an equal amount of charge on the island, no matter what the capacitance (and thus the 
distance) between them is. In other words, this result follows directly from using the "constant 
capacitance model" and neglecting screening. Although ref. M discusses the TLTS influence 
in terms of a general model where the potential rather than the charge of a TLTS is kept fixed, 
producing both a distance and a size dependence of the TLTSs, it is not in contradiction to 
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the present treatment. We rather consider the range of small charge displacement and small 
distance from the SET island here (top left part of fig. 56) in ref. ||). These conditions 
are reasonable and have already been concluded in Both methods yield corresponding 

results in case of short TLTS-island distance. 

The values obtained for 1/7 easily translate into a characteristic tunnelling distance 5 
within TLTSs, using eq. (|). With D = 70... 100 nm, depending on the device, we obtain 
5 = (3.94 ± 0.19) nm, taking into account eight experiments with geometrically different SETs 
measured at 10 mK. This value agrees well with other experiments |6|-||,[l(|. 

The value we extract for S thus exceeds the thickness of the SET tunnel barriers, which 
is typically 2 nm or less. Hence the monitored trap states cannot reside in a tunnel junction 
itself, but must be located in the surface oxide or the substrate ||. The measurements 
suggest the tunnelling processes themselves not to be field-assisted, since the probability of 
a TLTS event does not depend on V g . Under these circumstances J«4 nm is a rather large 
value. In the case of tunnelling between trap states, however, the larger tunnelling distance 
can be compensated by lower barrier height, which can in fact be well below leV and still 
warrant electric field independence of the tunnelling process. Thus, the height of the energy 
barrier between the two TLTS sites is limited towards low energies by the maximum poten- 
tial difference (w 80 meV) and towards high energies by the necessary barrier transparency 
(1=3 0.7 eV). The neglect of multiple electron transitions within one TLTS in our model might 
be another reason for the large value of 5. 

Information on TLTS switching events at elevated temperatures is hard to access due to 
the increasing overlap between the main peak of the p(AQi) distribution and its tail, as the 
temperature increases. Measurements performed at temperatures up to 200 mK did not reveal 
any significant variation of I/7. 

The determination of I/7 (and S) also provides an estimate of the density of states of the 
TLTSs. The data presented in fig. [I] contain 572 Coulomb oscillation periods, 177 of which 
belong to the tail of p(AV g ). While the total energy range scanned by V g is 572 e/C g sa 1.8 V, 
the potential shift experienced between the two states of a TLTS is only S/D times this value. 
Hence, the average density of states is 



Owing to lacking information on the energy level of the whole TLTS and its inner structure 
(which might allow, for instance, more than one TLTS transition), -Dtlts does not easily 
translate into a density of states of the traps accommodating the TLTS. 

Defects with internal degrees of freedom can switch between metastable configurations also 
due to their interaction with a thermal bath or by photon excitations. Low-frequency excess 
noise even in high-quality devices is significantly generated by dynamic fluctuations of such 
defects |p^| . While a single fluctuator, also known as random telegraph fluctuator (RTF), 
shows a Lorentzian spectrum, an ensemble of RTF results in the typical 1// noise jl7J. It 
is reasonable and tempting to assume that the TLTS effects investigated in our studies are 
of the same origin as RTF generated noise. Due to the metastable nature of the TLTS, the 
dynamics of a fluctuator crucially depend on its energy threshold separation. However, this 
is only vaguely known from our experiments via the observation of Coulomb peak pattern 
hysteresis |1. Furthermore, it is not clear how to extrapolate the data, where we can resolve 
TLTS events with a large hysteresis, towards a small energy interval of fcgT, i.e for fluctua- 
tors of small time constants contributing to non-zero frequency noise. Therefore, a decisive 
conclusion on consequences of our results for dynamic fluctuations and noise is too speculative 
at the moment. However, his subject is under further investigation. 
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Conclusion. We have investigated charge transfer in two-level tunnelling systems 
(TLTSs). The study is based on measurements performed on metallic SETs. The charge 
transfer can be described in terms of a Poisson process, where the influenced charge Qi re- 
places the usual time variable. Thus, the prominent tail of the probability distribution function 
p(AVg) is governed by a single scale parameter I/7. We can associate this parameter with a 
length scale 5 of the TLTS. We find 5 = (3.94 ± 0.19) nm for the samples of different geome- 
tries under consideration. From this result we conclude that the SET tunnel barriers are in 
our case free of TLTSs. The fluctuators reside in the oxide around the island electrode or in 
the substrate. The height of the energy barrier separating two TLTS states is between 80 meV 
and 0.7 eV and the density of states of TLTS is found to be 2.1 (mcV) -1 . 

* * * 

We are indebted to Sergey V. Lotkhov for the excellent sample manufacture, as well as to 
Mathias Wagner and Ulrich Zulicke for thorough reading of the manuscript. 
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Fig. 1 - Determination of 1 /'y from a fit to experimental data. Equation ([]) is used to fit the tail of 
the probability distribution function p(AQi). The main plot is enlarged to emphasize the fit of the 
tail, whereas the left and right inset show the complete p(AV g ) and p(AQi) histograms, respectively. 
The dashed lines in the plots represent the fit, and the arrow indicates the limit of the fit range. 




Fig. 2 - Background charges in two-level tunnelling systems, a) Schematic drawing of the SET's 
gate (f) and island (4) electrode, which are separated by the distance D. Because the island is 
quasi-isolated, the source-drain electrodes are not considered for electrostatic arguments. A TLTS 
is located between island and gate, where a charge can tunnel a distance S from one site (2) to 
another (3) or vice versa, b) The island charge Qi displays arbitrary jumps at definite values of V g . 
c) Transition to uniform V g steps at arbitrary Qi values. 



